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Calculations of the hadronic Bc*'-mesons production performed in the framework of 
the perturbative QCD taking into account O(o4) Feynmann diagrams are presented. A 
comparison of the exact calculations with those based on the fragmentation model of 
b — > Be*' + X shows the large discrepancy between them. The exact calculations of the 
B^'-mesons production cross-sections as the function of pr at the energy of the FNAL 
Tevatron (y/s =1.8 TeV) are given. The predicted ratio of the vector to the pseudoscalar 
state cross-sections is about R ~ 3 instead of R ~ 1.4 for the fragmentation model. 



I. INTRODUCTION 



The production of hadrons with two heavy constituent quarks of different flavours was the subject of the 
theoretical studies over two the last years. The calculations of the Be -meson production cross-section 
in e + e~ annihilation at the Z° pole [1-3] show a principal possibility of the ground state observation of 
the be system as well as its higher excitations at LEP-I energies. It is much more possible to observe 
Be -meson in experiments at the hadronic collider where we can expect [4-6] a large number of events 
at the available luminosity of the FNAL collider. 

Another aspect of interest to the processes of the heavy quark system production consists in the fact 
that it offers a possibility to understand the mechanism of the hadronization in the framework of the 
perturbative QCD with the minimal assumptions concerning wave functions of the bound states. 

Let us enumerate the main results of the studies of the be system production mechanisms in different 
reactions: 

i) £>c*'-mcsons production in e + e~ annihilation at high energy (M B / s ^ 1) can be described in the 

framework of the fragmentation model with the simple factorized expression 

da B (») 

—^-=°bb-D- b ^ B M{z), (1) 

where z = 2Eb c /\/s, and Dr B ( r )(z) is the fragmentation function of b — > B c [1-3]. 

The ratio of the vector _B*-meson production cross-section to that of the pseudoscalar B c is 

instead of the expected value R ~ 3 that one can obtain counting the quark states. 

The production of the P- level states in the fragmentation model is about by an order of magnitude 
suppressed with respect to that for 5-states [7] . 

The contribution of the c-quark fragmentation into the final cross-section is about two orders of 
magnitude lower comparing with that of the &-quark. 

ii) The production of B C (B*) -mesons in the 77 collisions is described by 20 Feynmann diagrams [8-10] 

that can be split into three gauge invariant groups (6+6+8 diagrams, correspondingly). The first 
two groups describe the process of the bb and cc pairs production with the subsequent fragmentation 
into Be -meson. The contribution of the c-quark fragmentation in 77 collisions is enhanced due to 
the factor (Q c /Qb) A an d can't be neglected contrary to the case of e + e~ annihilation. 

The contribution of the 6-quark fragmentation can be sufficiently well described by equation (1) 
and the accuracy of this description increases with the growth of px value. The contribution of the 
c-quark fragmentation is in a quite drastic contradiction with factorized expression (1). The main 
contribution into -B c (-B* )-meson production at any pt comes from the remaining 8 diagrams of the 
recombination type. This results in another value of R > 3 different from that in e + e~ annihilation. 



iii) The hadronic production of Be* -meson requires the calculation of the full set of 36 Feynmann 
diagrams of the order of 0(oq) including fragmentation type diagrams. The cross section of the 
B C (B*) can be then obtained by the convolution of the gluon-gluon cross-section with the gluon 
distribution functions in the initial hadrons. This leads to the dominance of the region of small 

\/i of the sub-process gg — ► Be + c in the cross-section of B c production. Just in the same way 
as in 77 collisions the recombination type diagrams dominate in gg production of B C (B*) [6] and, 
consequently, the R value is close to 3 in the hadronic production as well. 

It is obvious that the experimental cuts on px values rises up the role of large s in the reaction 
gg — ► Be +X. In such a case we should expect the increase of the fragmentation mechanism contribution 
and, consequently, the simplification of all calculations. The calculations at large pt values have been 
performed in [5,11] where the cross-sections of B c (B*)-meson production for the FNAL Tevatron are 

presented. The authors of the cited papers used as the fragmentation function for b — > Be the value 
obtained in e + e~ annihilation for the S and P levels. Then, convoluting this fragmentation function 
with the cross-section of the bb pair production they calculated the distribution of B C (B* )-mesons as the 
function of pr- The leading logarithmic correction on the final state gluon radiation was also taken into 
account. 

In this paper, considering the full set of Feynmann diagrams of the order of 0{a i s ) we will show that 
the fragmentation approach fails, because not all the diagrams, which can't be neglected at large pr, 
were taken into account and there was some overestimation of the real phase space. 

We will give a more detailed analysis of the B C (B*) production at large px for the energy of the FNAL 
Tevatron (y/s = 1.8 TeV). 



II. HADRONIC PRODUCTION OF B { * ] WITH LARGE P T 

It seems obvious at the naive level that the dominating mechanism of the Be -mesons production at 
large pr should be connected with the fragmentation of heavy 6-quarks. The fragmentation approach is 
based on the assumption that the factorization of the production on the parton (b) production with large 
energy and its subsequent fragmentation into different B c states is valid. The differential cross section 
da/dpT, for example, for pp collision will have the form 

-^(PP -> H(p T )x) j dxidx 2 dzf i/p (x 1 , n)f j/p (x 2 , ^) x 

X -T~( i j -» Kpt/z) + x) x D- b ^ H {z,n), (2) 

where D(z,fi) are fragmentation functions b — > H; H = B C ,B*..., da/dpx is bb pair production cross- 
section and fi/A{ x T^) is parton i density function in the hadron A. This approach was used in [5,11], 

where the differential distributions da/dpr for Be -meson and its P-wave excitation were calculated. 
Nevertheless, it remains unclear what is the region where we can use expressions (2) and what can be 
the role of other subprocesses. Let us clarify these points in the framework of the same approximation 
- the Born- level diagrams of the order of 0{oq) [6]. The total number of diagrams in the subprocess 
gg — ► B c + X is 36 and only minor part of them is connected with the 6-quark fragmentation. 

We have already showed [6] that the main contribution into Bc^-meson production is connected with 

the diagrams of the recombination type, where both initial gluons dissociate into a pair of heavy quarks 

(*) 

and two of these quarks recombine later into B c -meson. We observed the dominance of such diagrams 
up to the energies of gg-collisions of the order of 1 TeV. Moreover, at the energies V§ ~ 30 GcV, that 
give the main contribution into hadronic production of Be -meson it is rather meaningless to consider 
the fragmentation mechanism at all, because the condition M 2 /s <C 1 is not valid and the pre-asymptotic 
terms in the definition of D{z) are large. So, we can use the fragmentation approach at large values of 
y/I only, but at large energies the fragmentation contribution is small as compared with that of the whole 
set of diagrams. 

In Fig. 1 we present the results of the exact calculation for the cross-section of the subprocess gg — > 
B C (B*) + X as a function of c.m.s. energy in comparison with that calculated in the fragmentation 
approach - the cross-section of the bb pair production multiplied by the probability of the fragmentation 
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b — > B c and 6 — > _B* , 3. • 10~ 4 and 4.15 • 10 -4 , correspondingly. One can see, that in the regions where 
the condition s /M B 3> 1 is valid and expressions (2) can be used the contribution of the fragmentation 
is not dominant. Contrary, at small energies it gives overestimated values due to the fact that in the 
fragmentation calculation one uses two-particle phase space (2 — > 2) instead of true three-particle space 
(2 — > 3) in the partonic subprocess. 

The same conclusion can be drawn from Fig. 2, where we present the differential cross-sections da/dpr 
for B c - and B*-mesons in comparison with the fragmentation mechanism at the energy Vs = 100 GeV in 
gluon-gluon collision. In the case of £? c -meson production we observe the saturation of the fragmentation 
at the values of pt > 30 GeV, while in the case of -B*-meson this saturation is postponed to the very 
edge of the phase space (pr > 40 GeV) . 

To stress the problems with the fragmentation picture let us consider the differential cross-section 
da/dpr for B c - and i?*-mesons at the energy of gluon-gluon collision = 20 GeV (Fig. 3). One can 
see that incorrect calculation of the phase space leads to the fact that the cross-section calculated by 
fragmentation expressions (1) exceeds that obtained from the full set of diagrams. 

The conclusion one can draw from the above considerations is quite simple: at the currently reachable 
Pt the fragmentation approach does not work and all the diagrams contributing into Be mesons should 
be considered. 

It is important to underline another point. The full and the fragmentation calculations give different 
predictions for the ratio of the cross-sections of B c and B* productions. The full calculation gives the 
value of R ~ 3, while the fragmentation model predicts another value R ~ 1.4. 

In Fig. 4 we present the differential cross-section da/dpr for B c and -B*-mesons calculated from gg 
cross-section convoluting it with structure functions of the initial protons for the energy of the FNAL 
Tevatron (y/s = 1.8 TeV). We take the structure functions from [12] fixing virtuality at the value of Q <~ 10 
GeV, because we consider only those Born diagrams contributions, which approximately correspond to 
those virtualities (a s ~ 0.2). The running virtuality and the strong coupling constant can be used when 
one considers next-order radiative corrections to this process. 

Two curves in the same figure correspond to the contributions calculated from expression (2) using 
fragmentation functions Di^ Bc (z) and Di^ B ,{z) obtained at the same order of perturbative expansion. 
One can see that those curves do not coincide with the exact calculations in the whole interval of pr- At 
small values of px one can observe excess of fragmentation predictions, while at large values of px they 
are lower. 

Taking into account the experimental cuts on the pseudo-rapidity (\y\ < 1) one can calculate more 
realistic pr distribution (Fig. 5). Such a cut reduces the value of cross-section by a factor of 3. It is 
interesting to note that in the px interval considered the value R = ctb*/&b c is about 3. Taking into 
account all S'-wave excitations the cross-section of the i3 c -meson production IS(B C + B*) + 2S(B C + B*) 
with pt > 5 GeV is about 3.3 nb that gives for Run lb of Tevatron with the integrated luminosity 
100 -j- 150 pb _1 about 3.3 5.0 • 10 5 S c -mesons. 

III. DISCUSSION 

There are two different approaches to the calculation of the hadronic production of Bc^-mesons at 
present. The first one is based on the calculation of all the diagrams of the order of 0(a 4 ) in the 
perturbative QCD. There are three publications with such calculations, which give different predictions 
for cross-sections [4, 6, 13]. In our previous work [6], where we considered the hadronic B^ production, 

we had omitted the color factor l/\/3 in the wave function of Sc^-meson that increased our prediction 
by a factor of 3. Correcting this point and using the same constants as in [4] we have obtained a good 
agreement at the level of the gg cross-sections. 

The second approach is based on the usage of the fragmentation model [5,11], and, as we have demon- 
strated above, it gives incorrect description of the cross-section underestimating it at large energies of 
gg collisions and overestimating at small values of energies. So, it looks quite strange that the authors 
of [11] when using fragmentation approach and convoluting the resulting expressions with the structure 
functions have obtained the cross-sections, which are analogous to ours (see Fig. 5). 

Finally, we would like to note: 

1) In the case of the states with equal mass heavy quarks (say, ip production) the fragmentation approach 
c^> ip does not also describe the exact result obtained by the calculating of all the diagrams of the 
order of 0{a 4 s ). 
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2) When one quark becomes light (u, d, s instead of c-quark) in the region of large pt the recombination 
diagrams become dominant. 

All these points we will consider in our next publications. 
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FIGURE CAPTIONS 

Fig.l. Gluon cross-section in nb for the B c (white triangle) and B* (black triangle) production. The 
fragmentation model predictions for B c (dashed line) and B* (solid line) are given for comparison. 

Fig. 2. da/dpT distribution in nb/GeV of B C (B*) production (histograms) in gg collisions in comparison 
with the fragmentation approach (curves) at the energy of 100 GeV. 

Fig. 3. The same as in Fig. 2, but for the energy of 20 GeV. 

Fig. 4. Differential cross-section da/dpx for B C (B*) in pp-collisions at yfs — 1.8 TeV in comparison with 
the fragmentation approach. 

Fig. 5. The same as in Fig. 4, but with cut \y(Bi*^)\ < 1. The curves correspond to the predictions from 
[11]. 



4 



xlO 



S5 



0.2 - 



0.175 



0.15 



0.125 



0.1 - 



0.075 - 



0.05 - 



0.025 








Fig.2 



GeV 



xlO 




0.075 - 



0.05 - 



0.025 




Fig.2 



GeV 



